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Abstract 



Semiclassical approach has been developed for the one-dimensional interact- 
ing fermion systems. Starting from the incommensurate spin density wave 
(SDW) mean field state for the repulsive Hubbard model in ID, the non- 
Abelian bosonized Lagrangian describing the spin-charge separation is ob- 
tained. The Berry phase term is derived from the chiral anomaly, and we 
obtain the massless Tomonaga-Luttinger liquid in the single chain case while 
the spin gap opens in the double-chain system. This approach offers a new 
method to identify the strong-coupling fixed point, and its relation to the 
Abelian bosonization formalism is discussed on the spin gap state. The gen- 
eralization to higher dimensions is also discussed. 
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I. INTRODUCTION 



Correlation effects in one-dimension are characterized by the large quantum ffuctuations 
as well as the strong coupling nature of the interactions. As for the spin chains, it has 
been revealed that the quantal phase, i.e.. Berry phase, plays an essential role to determine 
the structure of the low energy spectrum Away from half-filling, the system becomes 
a Tomonaga-Luttinger liquid (TL-liquid) whose excitaions are exhausted by the collective 
modes, i.e., bosonic degrees of freedom 0. This TL-liquid is conveniently described by the 
bosonization method where the Hamiltonian is decoupled into spin and charge parts. All 
the above features are considered to be beyond the scope of the conventional mean field -|- 
RPA theory which gives invalid results even qualitatively in ID. Hence the interacting ID 
Fermion systems are regarded as a laboratory to study the correlation effects and test many 
theoretical techniques including the field theoretical methods and renormalization group 
(RG) [0. 

Another aspect of current interest is the spin gap formation and the possible supercon- 
ductivity in the double-chain system [§,0. This system can be regarded as the ID realization 
of the short range resonating valence bond (RVB) scenario proposed for the high-Tc super- 
conductors p. For the half-filled case, i.e., without the charge degrees of freedom, several 
authors have studied the double chain systems [|§-§[, and it is concluded that the spin gap 
opens both for antiferromagnetic (AF) and ferromagnetic (F) interchain exchange couplings. 
The mechanism for the spin gap is essentially the same as that proposed by Haldane for 
S = 1 AF Heisenberg chain 0. As for the doped case, mean field type theory [^, numerical 
works ||To|JTT|| , and analytical studies [|T^-[l^ have been done and the persistence of the spin 
gap even away from the half-filling has been claimed. It is noted that the fixed point for the 
spin gap state is outside the reach of the perturbative RG, and some assumption about the 
nature of this strong-coupling fixed point is nescessary to study its physical properties. 

Lastly the generalization of bosonization and Tomonaga-Luttinger liquid to higher di- 
mensions has been discussed intensively |T^, and it is highly desirable to develope a theoret- 
ical framework which is not restricted to one-dimension but correctly reproduce the results 
in one-dimension. It would clarify the role of dimensionality, nesting condition etc. in the 
physics of strongly correlated systems. 

In this paper we develop a semiclassical approach to the interacting ID fermions. This 
approach is reliable in higher dimensions, but has been believed to be a poor approximation 
even qualitatively in low dimensionality, especially in one-dimension. However we will show 
in this paper that in spite of the large quantum fluctuations the mean field theory is a good 
starting point even in ID providing a clear physical interpretation of massless Tomonaga- 
Luttinger liquid as well as the strong-coupling fixed point with the spin gap. 

For the model with repulsive interactions the appropriate mean field solution is the spin 
density wave (SDW) state |TB|. We identify the low lying collective modes, i.e., Goldstone 
modes, around the mean field solution and derive the effective Lagrangian for these modes. 
By introducing the rotating coordinate system of the spin [T7,TS], one can study even the 
case where the amplitudes of the collective modes are large and the long range ordering is 
absent. In this formulation, the SU{2) connection and hence the gauge field is naturally 
derived. This gauge field is interacting with the ID fermions and the chiral anomaly occurs, 
which results in the Berry phase term. The charge part, on the other hand, is described by 
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the U{1) phase which originates from the phason degrees of freedom [0, and is decoupled 
with the gauge field mentioned above. Then the effective Lagrangian is nothing but that 



of the non-Abelian bosonization with spin-charge separation |1T9| . In terms of this analysis, 



we obtain the massless Tomonaga-Luttinger liquid for the single-chain case while the spin 
gap opens for the double-chain case even away from the half-filling. This spin gap state is 
identified in terms of the Abelian bosonization. By analyzing the correlation functions, it 
is found that the low energy dynamics is described in terms of the bipolaron model pO 



These results in one-dimensional systems are based upon the collective-mode description of 
the low energy physics. This is possible only when all the individual excitations have a gap 
induced by the order parameter of the starting long range ordering, e.g., spin density wave. 
In higher dimensions this is not the case in general. When the nesting condition is well 
satisfied, however, the effective action remains similar even in higher dimensions. When the 
nesting condition is not satisfied and some parts of the Fermi surface remain unbroken, the 
fermionic field enters into the effective action. This fermionic field interact with the gauge 



field, and the action becomes that of the slave-fermion scheme |18 



The plan of this paper is the foUowings. In section II the semiclassical approach with the 
rotating frame is developed for one-dimensional systems. In section III the strong-coupling 
fixed point with the spin gap is identified and analyzed in terms of the Abelian bosonization. 
A short version of this part has been already published by one of the present authors (N.N.) 
20| , and we add the comparison with the semiclassical approach. Discussions including the 



generalization to higher dimensions and some conclusions are given in section IV. 



II. SEMICLASSICAL APPROACH WITH ROTATING FRAME 

We start with the Hubbard model in ID. 

H = -tJ2{ctci+ia+h.c) - f^oJ2 cf^Cia + UY^ Tli^nn (1) 

with the standard notations. The Hubbard interaction is rewritten as 

Uui^nn = ^(riii + nil) - ^(cL'^a/JCi/?)^ (2) 

where ffa/s is the (a/?) component of the Pauli matrix a = {a^, Cy, az). Therefore eq.(2) means 
that the repulsive interaction favors the formation of local spin moment Si = ^cl^aaisCijS. 
By introducing the vector Stratonovich- Hubbard field (^(r*, r), the partition function of the 
system is given by 



Dc^DcD^ exp 



/3 

- / Ldr 







(3) 



where 



1,(7 
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The mean field state corresponds to the static saddle point for the <^i(r) integration, i.e., 



U 



(s.p.) 



2me^ cos{2kpXi 



(5) 



where kp is the Fermi wavenumber and the spin polarization is assumed to be along the 
z direction, and the amplitude m is determined by the self-consistency condition. Because 
the amplitude fluctuation is massive, we neglect it and consider only the variation of phase 
6 and direction n {\n\ = 1) of the order parameter corresponding to the Goldstone modes of 
this SDW state 
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^(fiiir) = 2mni{T) cos{9{xi, r) + 2kFXi 



(6) 



where ni(r) and 6{xi, r) are assumed to be slowly varying, and we now employ the continuum 
approximation. 



dx 



m{n ■ c?)e~*^. 



m[n ■ a)e 
dr + ieAo + ivpidx + ieA^) 





' R' 




L 



(7) 



Here R 



R^ 



- ) represents the right going (left going) branch of Fermion with 

the linearized dispersion -^vpk measured from the Fermi wavenumber -^kp^ and the external 
electromagnetic field {Aq,Ax) is introduced. We introduce the polar coordinate (C,''?) as 
n = (cos sin sin sin cos C)- Correspondingly we employ the rotating frame whose pr- 
axis coincides with n [^,0. Explicitly this rotation can be accomplished by introducing 
the SU{2) rotation matrix 



U 



(8) 



where 



e*(^-'7)/2 cos§ 



3j(fe+r;)/2 



COS 



(9) 



is the spin 1/2 spinor and h is corresponding to the gauge choice. By using this ?7, one can 
easily obtain Ua^U'^ = n-a. Then the cross term in eq.(7) can be written as mR^{n-a)e^^L = 
mwazL where R = e~^^^'^U^R and L = e'^^^f/^L represents the fermions in the rotated 
frame. The phason field 6{f, r) appears as the chiral gauge transformation. Hence we 
expect the chiral anomaly. 

The chiral anomaly appears as the result of the gauge invariant regularization of the 
diverging integrals. One useful way to treat this anomaly is the Fujikawa's method |^ ( 
see Appendix). By using this method and also taking care of the spin degeneracy factor of 
2, the Jacobian for the change of integral variables is obtained as 



(10) 
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where E{x, r) = drA^ — dxAo is the electric field. The above equation describes the accel- 



eration of the phason in terms of the electric field ||22[. Besides this Jacobian factor, the 
Lagrangian in eq.(7) becomes 



dx 



9_ - ^d^e + ieA^ + Ud-W, 



maz 



d+ + ^d+e + ieA+ + Ud+U^ 





' R 




L 
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where d± = iivpdx and A± = AQ±ivFA^. The explicit expression for Ud±W is given as 



Ud±U'^ = -{—d-i-rism(cosb + d±(smb)a^ 
+ -[d±rj sin (siiib + d±(cosb)ay 
+ -{d±r]cosC + d±b)a^ 



(12) 



where the first and second terms represent the spin wave degrees of freedom. 

The coefficient of az, on the other hand, describes the connection between the neighboring 

rotating frame. Hence we regard it as the U{1) gauge field, i.e.. 



a± = ao ± ivpttx = d±rj cos C + d±b. 



(13) 



The "electric field" e is related to the Skyrmion (instanton) number density n ■ {drU x d^n) 
as e = drdx — dxCLQ = —n ■ [drU x d^n). From eqs.(ll) and (12) there are two massive Dirac 
fermions coupled with the internal gauge field a^, i.e., up spin electron with mass m couples 



with +|a^ while down spin electron with mass —m couples with 
action S[a^ for is given as 



Z{{aj2},m)Z{{-a,/2], 



-m] 



Hence the effective 



(14) 



Because we are employing the Euclidean formalism and calculate the partition function, the 
complex phase factor is coming from the contour integral / A{fi)-dri in the first quantization 
formalism. Then the relation Z({— a^/2}, m) = [Z({a^/2}, m)]* holds. Now the chiral 
anomaly is again relevant. By the chiral gauge transformation R e'^/^-R, L — > e^'^^^Z with 
^ = TT, the sign of the mass is inverted, i.e., m —m. The Jacobian for this transformation 
is J = exp[— I / dxdre^x, r)] = explinQ], where Q = J dxdrn ■ {don x dxn)/4:7i is an integer 
called the Skyrmion number. This relation together with eq.(14) gives 



(15) 



which states that the negative interference occurs between the topological sectors with even 
and odd Skyrmion numbers. This topological term was first derived by Haldane for the 
Heisenberg antiferromagnet summing up the Berry phase terms of the individual spins 0. 
Our derivation generalizes his results to (i) itinerant antiferromagnet and (ii) incommensu- 
rate case. We have shown that the chiral anomaly for describes the acceleration by the 
electric field E while that for the topological term. One may wonder the validity of the 
continuum approximation employed to derive eq.(15). To check this point we have studied 
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the Berry phase for the adiabatic change of the gauge field for several noninteracting 
electron systems in ID on a lattice. We found that, although the phase of Z({a^/2},m) 
depends on the model, the relations Z{{a^/2},m) = exp (i7rQ)Z({a^/2}, — m) and eq.(15) 
are always valid. That is, while the Berry phase itself is regularization-dependent, our 
calculations leading to eq.(15) is universal. 

We now proceed to complete our derivation of the effective Lagrangian for ip and n. For 
that purpose, we expand the effective action in terms of the derivative [^. We employ 
the gauge invariant regularization scheme and expand with respect to d±6 and the spin 
wave part of Ud±W neglecting the derivative terms of A± and a±. This procedure can be 
performed by expanding the Trln with respect to these quantities, and the lowest order 
term is given by 



J2Tr[Goik,uj)Vf 



(16) 



where 



Goik,uj) 



+ (vpky + 



iuj + vpk, mcTz 



ma. 



ILO 



Vpk 



(17) 



and 



V 



0, 





\d+e + Ud+U^ 



(18) 



Becuase V already contains the derivative, we regard V as if it is a constant matrix to obtain 
the second order terms with respect to the derivative. Another comment on the integral in 
eq.(16) is that we preserve the Lorentz invariance, i.e., the symmetry between uj and vpk. 
Then the final result is 



S 



dxdr 



-i—OE 
27r 



1 



J Attvf 



(19) 



which is the sum of the standard phason Lagrangian plus the nonlinear a model with topo- 
logical term. We have considered the slowly varying 6 and n without assuming the small 
amplitude of these fluctuations. 

Some remarks are now in order. In our Lagrangian eq.(19) there is no U appearing, which 
means that eq.(19) describes the noninteracting fermions. This drawback can be easily fixed 
by taking into account the forward scattering, i.e., g2 and terms in the terminology of 
^f-ology, as follows. Note that the density of the right-going (left-going) fermions pr (pi) is 
given by pn = d^O / {Ativf) {pl = d-0 / {Ativf))- The forward scattering term is U{pR + pif 
which is translated to U + d-OY / {AtivfY, which should be added to eq.(19). Therefore 
we obtain the action related to 9 as 



Sft 



drdx- 



Attv 



F I 



U 



2ttvf 



(20) 
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This action is rewitten in terms of the charge exponent Kp and the renormahzed charge 
velocity Vp as 



So = drdx 



A-nVpKp 



(21) 



where Kp = {1 + U/27ivf)~^^'^ and f p = vpil + U /27tvf)~^^'^ ■ Thus the Hubbard interaction 
U manifests itself through the forward scatterings in the renormalization of the charge 
exponent Kp and charge velocity Vp. The spin part, on the other hand, is not affected by U. 
These features are in consistent with the exact analysis at least in the small U limit where 
our approach is expected to be valid. 

Now the Lagrangian is actually equivalent to that of non-Abelian bosonization which 
respects the U{1) x SU{2) symmetry [0. Note that the gauge field representing the 



Skyrmion number density is not coupled with the charge phase 9. The actions for spin and 
charge are completely decoupled. In the case of attractive U Hubbard model, the appropriate 
starting mean field solution is that of the singlet superconductivity where the gap is opened. 
Therefore only the charge degrees of freedom is left which is described by the Josephson 
phase which is canonical conjugate to the phason field described above and the Lagrangian 
of essentially the same form as in the repulsive case is obtained. The charge exponent Kp, 
however, is larger than unity with the dominant superconducting fluctuations. 

The periodicity of the lattice has been neglected thus far. This corresponds to the 
commensurability pinning of the phason. At half-filling it is well known that the phason 
degrees of freedom is absent because Q = 2kF is equivalent to —Q = — 2fci?. Therefore the 
system is insulating, and only the spin degrees of freedom, i.e., n are left as the low lying 
modes. This is exactly the Mott insulator in one- dimension. For other commensurability, 
e.g., quarter filling case, the commensurability pinning can be described in the following way. 
Let the Q = 2kF = G/N with G being the reciprocal lattice vector and an integer. Because 
the charge density induced by the SDW is due to the second order effect and proportional to 
cos(2Qa; + 26), 2MQ = {2M/N)G {M,N: integers) should be some multiple of G in order 
that 2M-th order energy contributes. In this case there appears a commensurability energy 
as 

j dxV2M cos[2Me{x) + a], (22) 

where V2M is some potential energy and a is some phase. The integer 2M is no less than 4, 
and V2M is irrelevant when the exponent Kp is near unity, i.e., weak coupling, in which we 
are interested in this paper. 

We now apply our semiclassical method to the double-chain (two-band) model. Our 
interest here is the possible coexistence of the metallic conduction and the spin gap in the 
models with repulsive interactions only Now consider the double-chain Hubbard model 
coupled with the interchain hopping t±_ and the exchange interaction J. First consider the 
case of t_L = 0. We can follow the discussion above and the saddle point configuration (^(^-p-) 
is ip'j^^ = ±<^P'' oc miz cos{2kFXi) depending on the sign of the exchange J. Here A and 
B are chain indices. The effective action for the phason and the spin wave is the same as 
eq.(19) for each chain and there appears the interchain coupling. 
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= -J j dxn^^^ (x) ■ n(^) (x) cos(^(^) (x) - 0^^) (x) ) . (23) 

Because of this coupling the low lying collective modes also satisfy this relation = ^^^'^ = 
, i.e., n = n^"^^ = ±n^^^ and 9 = 9^^^ = 9^^^ and the effective action for the phason and 
the spin wave is again the same as eq.(19) except the change of the Berry phase term, i.e., 
iTcQ becomes 2i7rQ ( for J > ) or ( for J < ). In both cases the spin gap opens and the 
phason mode gives the metallic conduction. Because of the spin gap, this state with the spin 
gap is stable against at least the small interchain hopping t±. Considering that the energy 
gain due to the spin gap formation is of the order of \J±\ while the kinetic energy of the 
doped carriers ( concentration 5 ) is of the order of 6t or 6t±, the spin gap state is expected 
to be stable as long as 6 max{t,t±) << \J±\. In the next section we study the various 
correlation functions in this spin gap state. For this purpose, the semiclassical method 
presented above is not convenient because the canonical conjugate field 9- to the phason 
field 9 has been already integrated over when the effective action is derived. However the 
considerations above uniquely determine the strong coupling fixed point which is translated 
into the language of Abelian bosonization. The effect of the interchain hopping tj_ will be 
again discussed and it will be shown that it is irrelevant at the strong-coupling fixed point 
corresponding to the spin gap state. 



III. SPIN GAP STATE IN ABELIAN BOSONIZATION 



In this section we relate the spin gap state discussed in the previous section to the 
strong-coupling fixed point in the Abelian bosonization scheme. We start with the following 
Hamiltonian in the Abelian bosonization formalism IHI. 



H = Ha + Hb + Ht. + Hj 



where 

Hi = Vn 



dx 



4:TTrjp \ dx 



d9 



Wx2 



+ I dx 



4:-Krjp \ dx 



(i)2 



with i = A,B being the chain index and 

Ht^ = -tL J dx\ij^^J\x)i^^£{^) + A^J\x)^Z\x) + h.c. 



(24) 



(25) 



(26) 



— J± dxSAix) ■ Ssix) 



(27) 



9^\(f)^J^) is the phase variable describing the charge (spin) degrees of freeedom, and Pj!''(M^*^) 
is its canonical conjugate momentum. P+''(M+'*) is related to the phase variable 9^l\(j)^l'^) as 



P 



de 



-/27r (M 



ii) 



-^1 \±y±j^ — -—^/2tx). In terms of these phase variables, the field operators 
of the electrons and the spin are represented as 



8 



exp 



(28) 



and 



exp 



= — e-'*-^[cos0? + cos(^;!^ + 2kFx)], 



ira 



= ^e'<^"[cos0? + cos(^? + 2kFx)], 



(29) 

(30) 
(31) 



(32) 



2tt vra 

Using these expressions the interchain interactions Hj^ is written in terms of 6^) and (f)^^ as 



7ra 



(ix<l cos((/)^'' 



cos^V^^ cos^f ^ + - cos(^i^^ - e?^' 



'+ 



dx 



(vra)' 



cos 



cos < 



COS 







(33) 



where the rapidly oscillating parts with the wavenumber ^Akp does not contribute to the 
integral. 

It is obvious from eqs.(6) and (30)-(32) that the phason field 6 in the previous section 
is nothing but 6+ for each chain. The director of the spin n is related to the spin phases 
0±. Now we try to identify the spin gap state discussed in the previous section in terms 
of Abelian bosonization. The spin gap state is characterized by the fact that the relative 
fluctuations of the spin fields n's and charge ^'s on the two chains are massive, i.e., fixed. 

Let us first consider the case of half-filling with the charge degrees of freedom being 
quenched. This corresponds to fixing 6*^'' in eq.(33), and the system is described by only the 
spin phases 0's. Several authors have studied this problem and it is concluded the 

spin gap opens even for infinitestimal \J±\- This can be understood as the interchain singlet 
formation for antiferromagnetic (AF) J± and as the Haldane gap |^ for ferromagnetic (F) 
J±. In terms of the Abelian bosonization method these spin gap states are described as the 



massive phase of 

q(A) q{B) MA) q{B) 



{A) 



l{B) 



for both F and AF J_i 



PJ^. This is because 



l{A) 



+ 



l{B) 



and cosff/)^^ 



) 

(B) 



+ (p)^' and 

gives rise to cos(0^^ — 0^'*) and S^-^^S^-^^ gives rise to cos 
. It is impossible to fix both the canonical conjugate pair 
and 0+"^ — 0\"^ and the spin gap state is realized by fixing — 0^'' and 

From the above construction, the relative direction of S^^^ and 5*^^^ is fixed, i.e., < 
S!/^^Sli^^ >, < S^-^^S^^^ >, and < S^-^^S^-^^ > are all nonzero, while there appears no spin 

moment ( < S^^^ >= (< S^^'^ >= 0). Therefore this state is identified as the spin gap state 
for n^"^) and n^^^ in the previous section where the relative phase of them is fixed. 
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Now let us turn to the doped case. The first question in this case is whether the spin gap 
survives the doping or not. When the interchain interactions are treated perturbatively, the 
simple power counting arguments give the following conclusions. As can be seen in eq.(33), 
the cosine type interactions in Hj^ discussed above are multiplied by cos(6'+^^ — 6'^-') and 
their exponents are increased by 2rjp. Therefore the interchain hopping Ht^ is more relevant 
than Hjx, and the conventional treatment is first to diagnalize Hti_ to obtain the bonding 
and anti-bonding bands and later to take into account the interactions. However we are 
interested in the situation where J±_ is reasonably large, and the electron number is near 
the half- filling, i.e., near the Mott insulator. Let 6 be the concentration measured from 
half-filling. Then the characteristic energy due to the phase fluctuation O'^^ is raax.{t^t±)8 
where t is the intrachain hopping. On the other hand, the stabilization energy due to the 
spin gap formation is of the order of |J_l|. Hence if max(t, tx)5 ^ |Jx|, it is allowed to 
consider first the exchange interaction Hj^ and later treat the interchain hopping Ht^ as a 
perturbation. 

In the case of t±_ = 0, the effect of the doping is summarized in the factor cos(^^^ — 0^^) 
which appears in eq.(33). If this factor gives the finite expectation value, i.e., the combination 
6^'' — 6^^^ is fixed and massive, the dynamics of the spin phases 0's remains essentially the 
same and the spin gap persists. On the other hand, if the spin gap persists, i.e., (f)L^ — 
and f^^'' + 0^^ remain fixed, cos(6'^'' — ^^^) alone is the relevant perturbation and 6'^'' —6'^^ 
is fixed, which means that the spin gap state is self-consistent. This state is identical with 
the spin-gap fixed point found by Khveshchenko and Rice |T^. This is also identified with 



the doped spin gap state in section HI where 9^^^ — 9^^^ is fixed and < n^"^^ -n^^^ > is nonzero. 
Because V9^+^ (V^f ^) is the slowly varing part of the charge density on the chain A{B) M, 
this means that the charge balance between the two chains is kept. If the carriers are doped 
with different concentration, 9^^ — 9^^^ is proportional to the cordinate x along the chain 
and cos(6'^'' — 6*^'*) does not contribute to the integral. The spin gap disappears immediately 
and it is expected that the system is described as the massless Tomonaga-Luttinger liquid 
in this case . 

This spin gap state has stability against small interchain hopping t±_ because Ht^ is 
the irrelevant perturbation in this case, i.e., its expectation value vanishes and correlation 
function decays exponentially. Also it should be noted that the two chains need not be 
exactly equivalent. If the chemical potential difference A/i between the two chains are small 
enough compared with | the stabilization energy of the gap formation prefers the equal 
concentration of the holes on each chain. 

Now we consider the physical properties of the doped spin gap state. As has been 
discussed above, this state is characterized as the fixed and massive state of the combinations 
0+^^ + (t>+ \ — and 9^_^'' — 6*^^ In this state the canonical conjugate fields to those 
fixed ones, i.e., 9^f -9^^\ 9i^H9^^\ and 9^^^ -9^^\ are uncertain and their exponentiated 
operators have zero expectation value and exponentially decaying correlation functions. The 
only massless mode is then + 9^^\ i.e., the in-phase charge fluctuations. From these 
considerations the following conclusions are derived immediately . 



(1) Among numerous order parameters constructed as the product of the two fermion fields , 
only the interchain singlet superconductivity Ogg given below shows the power-law behavior 
while all the others show exponential decay. 
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(34) 

(2) When one goes further to the product of four fermion fields, the 4kp charge density wave 
O^^,^^ given below shows the power-law behavior. 

qABAB _ , (A)t , (i3)t , (A) , (B) 
1 \2 



2Tia 



exp 



(3) The exponents Ksc and Kcdw for the two order parameters defined in eqs.(34) and (35) 
respectively are determined by the only massless field + ^ , and satisfies the duality 
relation Ksc ■ Kcdw = 1 HH- And it is noted that the wavenumber Akp is corresponding 
to the "2/cir" of the spinless hard core boson. Note that the 2kF charge density wave of the 
fermion does not show power-law correlation. Therefore this spin gap state in the double 
chain system is distinct from the Luther-Emery type state in the single-chain state P7| , p5| . 
The Luther-Emery state for negative U Hubbard model is understood as the fluctuating 
singlet superconducting state as has been discussed in sention II. These two facts means 
that the low energy dynamics of the spin gap state is described as the bipolaron along the 
lung of the ladder. This is also consistent with the above conclusion that the spin gap 
persists only when the two chains are doped equally. Recent numerical study of the double- 



chain Hubbard model |jrT[ suggests that Ksc is 2 which means Kcdw = 1/2, i.e., the CDW 
is dominating. However we have more chance to obtain superconducting state in the t-J 
ladder model because the naive expectation is that the product of two fermion operators 
have small scaling dimensions than that of four fermion operators. In this picture the charge 
carriers in the spin gap state is the charge 2e bipolaron. 



IV. DISCUSSION 

In this section we discuss the implications of the above results for higher dimensions. 
The Stratonovich-Hubbard transformation ( eqs. (l)-(4) ) and the rotating frame method ( 
eqs.(8)-(9) ) can be used in any dimensions. In eq.(4) there occurs a competition between the 
kinetic energy t and the local magnetic field 2U(fi/3. The former prefers the nonmagnetic 
state due to the Fermi degeneracy, while the latter induces the magnetic moments. The 
criterion for the weak and strong correlation is the relative magnitudes of these two, i.e., 
the bandwidth W ^ zt {z: number of the nearest neighbor sites) and the average strength 
of the local magnetic field F = 2U < \(pi\ > /3. 

When F << W the laboratory frame with the fixed spin axis is the appropriate coordi- 
nate in which the kinetic energy Hamiltonan is easily diagonalized. The local magnetic field 
is a small perturbation which slightly induces the magnetic moments. This can be treated 
as a weak scattering of the fermions due to the local magnetic field. The density of states 
remains essentially the same as the noninteracting case. 

When F » W, on the other hand, the fermion spin at site i is forced to be parallel 
to the direction of the local magnetic field Hi = ipi/\0i\. Therefore one must take into 
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account the local magnetic field first by employing the rotating frame whose 5'^-axis coincides 
with Hi = (cos?7 sin sinr] sin cosC). This can be done similarly to eqs.(8) and (9) by 
introducing the 2x2 unitary matrix Ui. 



JJ. ^ g-%o-z/2g-jCiO-y/2gi6,(T^/2 



Zil, 



(36) 



where 







-Zii. 





2 J 



Then we define the fermion field in the rotated frame as 



UiC = Ui 



(37) 



(38) 



The local "magnetic field" is always along the +z direction for the new fermion c. In this 
case the density of states for the fermions and c is split into two according to their spins 
with the gap of the order of F. This is nothing but the upper and lower Hubbard bands 
with the bandwidth of the order of W. Note that we do not assume the magnetic ordering. 
For each Hubbard band, the number of available states is that of the lattice sites N. Then 
for the half-filled case, only the up-spin band, i.e., lower Hubbard band, is occupied, and 
the system becomes the Mott insulator. When the holes are doped, there appears the small 
hole pocket near the top of the lower Hubbard band. In this case the upper Hubbard band 
is irrelevant for the low energy behavior of the system, and we neglect it for the moment. 
Then the system is described by the spinless fermion /j = and its conjugate. The original 
electron field is given from eq.(38) as 



(39) 



which is exactly the slave- fermion decomposition [|T^] . The low energy dynamics is described 
in terms of the spinless fermion for charge degrees of freedom and the spin field n (or Za). 
In this case the gauge field is coupled with the spinless fermion. 

The above discussion did not include the detailed band structure and the shape of the 
Fermi surface in the momentum space. When the nesting condition is satisfied, it is possible 
that the gap opens up even when F « W ii one chose the spatial pattern of (pi appro- 
priately. The one- dimensional system discussed above is the typical example of this case, 
where the "Fermi surface" is the two points with perfect nesting. We concentrate on the 
Fourier components of ifi near = ±2kF, and the low lying excitations are exhausted by 
the collective modes. These collective modes are nothing but the variables describing the 
Tomonaga-Luttinger liquid, and the spin-charge separation occurs when the forward scat- 
terings are taken into account. Therefore it represents the effects of the correlation on the 
collective coordinates. When the nesting condition is well satisfied, this scenario remains 
basically the same even in higher dimensions. Of course the long range ordering gener- 
ally occurs in higher diemnsions and the collective modes can be treated in terms of the 
standard random phase approximation (RPA) when the fluctuations are small. However, 
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near or above the transtion tempcratTirc, or in the case where the long range ordering is 
suppressed by some reason, e.g., frustration in the interactions, the fluctuations are large 
and one had better describe the system in the rotating frame similarly to the case of strong 
correlation. Therefore the system with the nested Fermi surface share some features with 
the strongly correlated system when the long range ordering is absent, and will contribute to 
the understanding of the role of the collective modes in the physics of the strongly correlated 
systems. 
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APPENDIX A: FUJIKAWA'S METHOD FOR CHIRAL ANOMALY 



In this appendix we give a brief description of Fujikawa's method apphed to the 
Dirac fermions in (l+l)-dimensions. Here we take the unit vp = 1 and let r = (r, x) be 

-R{r) 
L{r) 

[L'l"(r), i?''"(r)]. We take the chiral representation where 75 = a^. Then the chiral gauge 
transformation is given by 



the two-dimensional coordinate. %l){r) 



is the two-component spinor, and il){r) 



(Al) 



We are interested in the Jacobian corresponding to the change of integration variables, i.e., 
chiral gauge transformation given above. For this purpose let us expand ip{r) in terms of 

0, D, 



the eigenfunction (pm{T^) of the Dirac operator D = 







with D± = d± + iA±. 



and 



(A2) 



(A3) 



The the functional integral Dil){r) can be replaced by Ildam- Then also the transformed ip' 
can be expanded in terms of ipm as 



m 

with being given by 

n n 

Coresponding to this change of integral variables 

Yijnda^ = [detCmn] n„(ia„. 
Assuming small a(r), the determinant is explicitly given by 

exp -i / (fra{r)^ipl^{r)'-i^ipri 
exp —i I (Pra{r)G{r) 



(A4) 



(A5) 



(A6) 



[detCmn] ^ = det 



(A7) 



The function G{r) is calculated as 
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G{r) ^J2'PU'>^)^5^n{r) 

n 



M— ►oo 

— lim Tr 



M— ►oo 



/ 



(A8) 



where we have introduced the gauge invariant convergence factor e-(-^/^)' to remove the 
ultraviolet divergence. Now we examine D'^. It is easy to obtain 



= {dr + lAof + {d, + lA^f + 75FC 



01 



(A9) 



where Fqi = drA^ — dxA-,- is the electric field. Then the function G{r) can be calculated as 
G{r) = ^rni^ Tr J exp [- ^ {k, + A.f/M'^ - F,,^,/M^' 



IJ,=T,X 



— lim 



2Ff 



01 



M^oo M2 J (27r)2 

— -^01 
2n 

Putting eq.(AlO) into eq.(A7) we obtain 



cxp 



Eik, + A,r/M' 



IJ,=T,X 



[detC„ 



exp 



Foi(r) 



2n 



(AlO) 



(All) 



Similar calculation can be done also for Dt/j and the same factor as eq.(All) is obtained. 
Then the Jacobian J for the chiral gauge transformation is given by 



J = exp 



Foi{r) 



lid rayr) 

. J TT 



(A12) 



which is Fujikawa's Jacobian. 
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